September 11, 1992 



LBL-32872 
UCB-PTH-92/33 



Meson Mass Splittings in the Nonrelativistic Model [] 



Richard F. Lebed 



Department of Physics 
University of California 
and 

Theoretical Physics Group 
Physics Division 
Lawrence Berkeley Laboratory 
1 Cyclotron Road 
Berkeley, California 94720 



Abstract 

Mass splittings between isodoublet meson pairs and between 0~ and 
1~ mesons of the same valence quark content are computed in a detailed 
nonrelativistic model. The field theoretic expressions for such splittings 
are shown to reduce to kinematic and Breit-Fermi terms in the nonrela- 
tivistic limit. Algebraic results thus obtained are applied to the specific 
case of the linear-plus-Coulomb potential, with resultant numbers com- 
pared to experiment. 



*This work was supported by the Director, Office of Energy Research, Office of High Energy 
and Nuclear Physics, Division of High Energy Physics of the U.S. Department of Energy under 
Contract DE-AC03-76SF00098. 



Disclaimer 



This document was prepared as an account of work sponsored by the United States Gov- 
ernment. Neither the United States Government nor any agency thereof, nor The Regents 
of the University of Cahfornia, nor any of their employees, maJjes any warranty, express or 
imphed, or assumes any legal liability or responsibility for the accuracy, completeness, or use- 
fulness of any information, apparatus, product, or process disclosed, or represents that its use 
would not infringe privately owned rights. Reference herein to any specific commercial prod- 
ucts process, or service by its trade name, trademark, manufacturer, or otherwise, does not 
necessarily constitute or imply its endorsement, recommendation, or favoring by the United 
States Government or any agency thereof, or The Regents of the University of California. The 
views and opinions of authors expressed herein do not necessarily state or reflect those of the 
United States Government or any agency thereof of The Regents of the University of California 
and shall not be used for advertising or product endorsement purposes. 



Lawrence Berkeley Laboratory is an equal opportunity employer. 



ii 



1 Introduction 



The splitting of the masses of mesons in an isospin doublet, sometimes called 
electromagnetic splitting, has traditionally been attributed primarily to explicit 
isospin breaking {i.e. rriu 7^ 171^) and differences between the charges of the 
valence quark-antiquark pairs {Qu 7^ Qd), with hyperfine, spin-orbit, and other 
effects neglected in comparison. Such a model serves to explain the observed 
splittings -K+ = 4.024 ± 0.032 MeV and D+ - D° = 4.77 ± 0.27 MeV, but 
has failed in light of the surprisingly small — = 0.1 ± 0.8 MeV. 

It is precisely this mass difference which has led to the proposal of a variety 
of models. Some of these [Q, |^, ^ ^ are based on the nonrelativistic model of 
hadron masses put forth by De Riijula, Georgi, and Glashow p soon after the 
development of QCD. Such models have the unfortunate tendency to predict 
numbers no smaller than B^ — B^ ~ 2 MeV, well outside the current exper- 
imental limits. Using more phenomenological models 0, one can obtain a 
smaller splitting in closer agreement with experiment. Nevertheless, it may seem 
odd that the usual nonrelativistic model, which works well for the D- and even 
the K- mesons, should fail in the case of the i?, which boasts an even heavier 
quark. 

The primary conclusion of this work is that it is possible to explain the mass 
splittings of heavy mesons {D and B, but not K) in an ordinary nonrelativistic 
model, as long as we take into account all corrections to consistent orders of 
magnitude, that expectation values of the mesonic wavefunctions in general 
have mass dependence, and that the running of the strong coupling constant is 
not negligible. 

In this spirit, the paper is organized as follows: in the second section we 
consider the problem of computing mesonic mass contributions in field theory. 
Then, in the third section, we demonstrate that the nonrelativistic limit of the 
field-theoretic result leads to kinematic terms and the Breit-Fermi interaction, 
exactly as stated in De Riijula et al. This is followed in Section 4 by an exhibi- 
tion of the full mass splitting relations for isodoublet 0~ and 1~ meson pairs, as 
well as (0", 1^) pairs with the same valence quarks. Section 5 discusses the ap- 
plication of quantum-mechanical theorems, including a very useful generalized 
virial theorem, to the problem of reducing the number of independent expecta- 



1 



tion values in the splitting formulae. These theorems are applied to the popular 
choice of a linear-plus- Coulomb potential in Section 6, with numerical results 
presented in Section 7. 



2 Mass Computation in Field Theory 

Typically, the computation of mesonic mass splittings in a nonrelativistic model 
is accomplished by starting with the Breit-Fermi interaction Sees. 38-42] 

Hbf = ^(aQiQj + kas) 

i>j 
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where rj, pj, mj, Sj, and Qi denote the coordinate, momentum, (constituent) 
mass, spin, and charge (in units of the protonic charge) of the zth quark, re- 
spectively; Tij = Ti — Tj] a and are the (running) QED and QCD coupling 
constants; and /;; = — 1(— |)isa color binding factor for mesons (baryons). This 
expression includes an annihilation term if qi = qj are in a relative j = 1 state. 
From this, one chooses the terms that are considered significant and then calcu- 
lates the appropriate quantum mechanical expectation values. We will pursue 
this course of action in the next section; however, this author feels that it would 
be worthwhile to first consider the derivation of this interaction for the mesonic 
system from the more fundamental field theories of QED and QCD, since this 
approach entails greater generality and may provide impetus for work beyond 
the scope of this paper. 

We first consider the question of the mass of a composite system from the 
point of view of the S'-matrix and interaction-picture perturbation theory. The 
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mass of a system, defined as the expectation value of the total Hamiltonian in 
the center-of-momentum frame of the constituents, receives contributions from 
both the noninteracting and interacting pieces of the Hamiltonian; the former 
gives rise to the masses and kinetic energies of the constituents, and the latter 
produces the interaction energy. Technically, the matrix element of the non- 
interacting piece in the interaction picture produces terms which contribute to 
interactions between renormalized constituents. Thus one may think of interac- 
tions between "dressed" constituents, a topic to which we return momentarily. 

Let us follow the method of Gupta to derive the interaction potential 
from the field-theoretical interaction Hamiltonian. We begin by writing the 
S'-matrix in the Cayley form 



S 



l + \iK' 



and expand the Hermitian operator K: 



(2) 



(3) 



The purpose of this expansion, rather than expanding S directly, is to preserve 
unitarity in each partial sum of S. The physical effect of this parametrization is 
to eliminate diagrams with real intermediate states from the S'-matrix expansion. 
Computing the terms i^^, one finds 



dtHUt), 



(4) 



where / indicates the interaction picture. Now observe that we may invent an 
effective Hamiltonian, H^jj, such that its first-order contribution is equivalent 
to the contribution from H-^^ to all orders. Thus, 



K 



dtHl^f(t). 



(5) 



The interaction energy is then 







^0 


if 





(6) 



with i^) and /^) actually the same state since the system is stable. 
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In our case, in which Hl^^ is composed of the interaction terms of QED 
and QCD, the lowest-order contribution is , corresponding to two interaction 
vertices: the exchange of one vector boson. It is easily shown that 



where the superscript (2) indicates second-order in the coupling constant, and M. 
is the usual invariant amplitude for the process. Eliminating the delta functions 
that arise in the right-most expressions we find 



Beyond second-order the relation between interaction energy and the invariant 
amplitude becomes less trivial, but nevertheless Gupta has shown that it can be 
done. However, we do not continue to fourth-order in this work, and henceforth 
suppress the (2) in the following. 

In general, Ai fi at any given order is represented by diagrams of the form 
indicated in Figure 1. The composite state is formed by superposition of the 
constituent particle wavefunctions in such a way that the desired overall quan- 
tum numbers for the composite state are obtained. For the mesonic system, 
A4fi is represented by the diagram in Figure 2, where the lowest-order interac- 
tion is the exchange of a single gauge boson. This class of diagrams allows for 
only the valence quark and antiquark (no sea qq pairs or glue), and thus would 
be a poor model if we chose these to be current quarks. Instead, the quarks in 
our diagrams will be constituent quarks, and the gauge couplings will assume 
their running values. In this way we can model the hadronic cloud, as well as 
renormalizations of the lines and vertices of our diagram, so that its particles are 
"dressed" in two senses. There is also an annihilation diagram if the quark and 
antiquark are of the same flavor. In this work we consider only the exchange 
diagram, since the mesons of greatest interest to us are those with one heavy 
and one light quark. 

The next step is to obtain the amplitude Aifi, in which the constituent legs 
are bound in the composite system, from the Feynman amplitude Ai [Figure 3] 
for the same interaction with free external constituent legs. To do this, we need 
only constrain the free external legs in a way which reflects the wavefunction 



(7) 
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and rotational properties of the meson state. In general, if the variables Zn are 
the degrees of freedom of the meson state |$), then we may write 

\^)=^dZn(t>{Zn)0{Zn)\^). (9) 

The function is an amplitude in the variables wavefunction; and O 

is a collection of Fock space operators which specifies the rotational properties 
of 1$). The integral-sum symbol indicates summation over both continuous and 
discrete Zn- In this notation, we obtain the result 

= ^ dzf dz, ^*izf)^{zi)f{zi, Zf)M{zi, Zf), (10) 

where f{zi,Zf) = (0| O"^ {zf)0{zi) |0) is a constraint function. We have written 
the energy contribution in this very general way in order to demonstrate the 
power of the technique. 

Now we apply this prescription to the usual case of Feynman rules. Then 
Zn are quark momenta, is the mesonic momentum-space wavefunction, and / 
specifies the spin of the meson, as we shall see below. The energy contribution is 
evaluated in the quark center of momentum frame {i.e. the meson rest frame), in 
which the relative momenta of the quark-antiquark pair, intitially and finally, are 
denoted by p and p', respectively. Fourier transformation of the wavefunctions 
from momentum-space to position-space yields 

AEcM = Jd^Xf Jd^Xi^*{xf)K{xf,Xi)tp{x,i}, 
where 

K(xf,Xi) = J d'^p' J d^p exp[i{p' ■ Xf — p- Xi)] f (spins) Ai{p' ,p, spins), 

spins 

and j d^xil)*{x)ip{x) = I. (11) 

As a technical point of fact, it is necessary to keep track of the normalization 
conventions used for wavefunctions, Fourier transforms, and Feynman rules in 
order to obtain the true convention-independent Ai?. As it stands, Eq. ^ locks 
us into a particular set of Feynman rule normalizations, which should be made 
clear in the following expression. The kinematic conventions are established in 
Figure 4. Then the Feynman amplitude for free external quark legs and a virtual 
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photon is 



M 



(27r)3/2 

VHXPi){-'iQnt.)vHfiPf) 




'121 



with Qqe^ replaced by for the gluon-mediated diagram. Note the use of 
hehcity rather than spin eigenstate spinors, which is done in order to implement 
a relativistic description of the mesons. In a nonrelativistic picture in which 
meson spin originates solely from the spin of the quarks (s-waves), spin-O(l) 
mesons have spin-space wavefunctions described by the usual singlet and triplet 
quark wavefunction Qq combinations: 



V2 ^ 



\ I spins. 



(13) 



The above expression remains true in a relativistic picture if we take the initial 
and final spin-quantization axes to coincide with the axes of relative momenta p 
and p', respectively, and then take T; i as helicity eigenstates. This is nothing 
more than the simplest nontrivial case of the Jacob- Wick formalism [1^. It is 
then a simple matter to write the constraint function for singlet (triplet) mesons: 



f{heUcities) = -^(^/.^i^^a ± Sh^iSn,^) ^(^/.^t'^h^I ± ^hfi^Hf^), 



(14) 



and so the object of interest is the constrained matrix element Aismg or Aitnp, 
which is the Feynman amplitude multiplied by the constraint function and 
summed over spins (or helicities). This is the object that is Fourier transformed 
in Eq. ^ 

In summary, mass contributions due to a binding interaction in a system 
of particles may be computed by writing down the Feynman amplitude induced 
by the interaction Hamiltonian, constraining the component particles to satisfy 
the symmetry properties of the system, and convolving with the appropriate 
system wavefunction. The specific implementation of this technique to spin-0 
and spin-1 mesons with constituent quarks in a relative i = state is described 



by Eqs. O, 12, and 14. 
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3 The Nonrelativistic Limit 



With the method for computing mass contributions in hand, we find ourselves 
with two possible courses of action. The first is to compute Msing or Mtrip in ^ 
fully relativistic manner, and then Fourier transform the result to obtain AEqm- 
The second is to immediately reduce the spinor bilinears via Pauli approximants, 
thus producing a nonrelativistic expansion. Let us explore both directions for 
the pseudoscalar case; the vector case is not much different. 

The relativistic result is noncovariant, because the energy contribution is 
evaluated specifically in the CM frame of the quarks. We see this refiected in the 
computation of the matrix element. For example, it is convenient to eliminate 
spinors from the calculation by means of relations like 

2^Uh{pA)Uh[PB) = I ^ — 1 ^ =, (15) 

h yj2mA[EA + mA) ^ yj2mB[EB + mB) 

and the explicit 70 is a signal of the noncovariancc. Once the spinor reductions 
and the resultant trace are performed, we find the expression 

Msing = -(Qge^ + glWr^, (16) 

where M results from the normalization factors, and T is the gamma-matrix 
trace. They are given by 

^ ^ ^^[i^.(i^. + M)i^Ki?; + M)e,(e. + m)e^(5^ + m)]-^/^ 
and 

r = 8 {{pi ■ Pi) [2efEf + 3{mEf + Msf + mM)] 

+ {pf ■ Pf) [2eiEi + 3{mEi + Msi + mM)] + {pi ■ Pi){pf ■ Pf) 

-{pi ■ Pf) [2EiEf + M{Ei + Ef + M)] 

-{Pi ■ Pf) [2e^ef + m{ei + ef + m)] + {pi- pf){Pi ■ Pf) 

-{pi ■ Pf) [mEi + Mef + mM] 

-{Pi ■ Pf) [mEf + Msi + mM] - (p, • Pf){Pi ■ pf) 

+ [-2mM{Ei-Ef){ei-ef) 
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+mM {m{Ei + Ef) + M{ei + £/) + mM) 
+2m^EiEf + 2M'^eief] } . 



(17) 



Also, 



{Pi-Pff = {ei-eff-{p~pr. (18) 

It is, in principle, possible to Fourier transform the product Aising of these 
unwieldy functions to obtain the full relativistic result for AEcm] this has not 
yet been performed. We can also perform the expansion of the energy factors 
in powers of ^, where all such momentum-over-mass quotients that occur are 
taken to be of the same order. 

However, this is unnecessary work, for if we require only a nonrelativistic 
expansion, there is a much faster way, namely expansion of the spinor bilinears 
via the Pauli approximants 

u(p')7m(p) 



/ /I {P + P') , X (P' -P) I \ , 
(X I + « 7, \X) + o 



2m 



2m 



/ /M , {p + p'f , (p' X 
(X I 1 + + ^ — \X) + o 



im^ 



Am^ 



(19) 



Using these expansions in Eq. ^ and taking \x) , \x') i^i helicity basis, we quickly 
find 



M 



sing 



{2Tif{p-p'Y 

41 ~1 



1 + 



{p + p'Y {p-p'T 



AmM 



m^ 



A 1 

+ 



mM M2 



+0 



(20) 



The gluon diagram has the additional physical constraint that the initial 
and final qq pairs are combined into a color singlet; this introduces an additional 
factor of — |. Then Fourier transformation of this result produces 



CM,sing 



A 
o 



r/^ 2mM 



-77 ^ rT + 



vr / 1 

2 Vm2 mM ' MV 



(21) 
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In comparison, the energy contribution from the Breit-Fermi interaction (Eq. 0) 
for a quark-antiquark pair of masses m, M in the CM reduces to 



(Hbf) = [oiQq - '^a)j 

^(^(p^ + r-(f-p)p)) 



r / ^ 2mM 



+ 



+ 



M2 mM 



q flavors ) 



-(- 

2 \r' 



^ I ^ 2mM/ 



[22) 



where G = ^ (^Sq ■ sq^, which is —1 for S = 0(1). Also, S = Sq + sq, and 5* 
is the AL = 2 tensor operator 



5*12 = 3(^1 • f)(a2 ■f)-ai- ^2- 



(23) 



For mesons with differently-flavored quarks in a relative £ = state, many of 
the terms drop out. Let us define 

C = <^-(p^ + f ■ (f ■ p)p)^ , 

D ^ (6'{fy 



(24) 



Then Eq. 52] becomes 



{Hbf) 



4 

aQq - -as 



B + 



1 TT / 1 J_ AQ \ 

2mM 2 Vm2 ^ ^ mM J . 



(25) 



and this is exactly Eq. ^ where Q = —1. 

We have been up to now considering only the contributions to the mass 
originating from the binding interaction due to one-gluon and one-photon ex- 
changes; there are, of course, also contributions from the kinetic energy of the 
quarks. Were we calculating these quantities in a relativistic theory, we would 
simply compute KE = (^^/m? -\- p^^. The square root may be formally expanded 
in norelativistic quantum mechanics as well, resulting in an alternating series 
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in (p^"). However, for large enough n in NRQM, these expectation values tend 
to diverge. For example, in the hydrogen atom, divergence occurs for s-waves 
at n = 3. Furthermore, if the system is not highly nonrelativistic, the inclusion 
of the {p^) may cause us to grossly underestimate the true value of the kinetic 
energy. The problem is that there is no positive (p^) term to balance the large 
negative (p^) term. For these reasons, we incorporate the alternating nature of 
the series in a computationally simple way by making the ansatz 

KE = ^m^+ (p2). (26) 

In order to evaluate the expectation values in the above equations, we will 
need to choose a potential. In the meantime, let us simply denote it with U{r). 
Then at last we have the mass formula: 

M„ = + (P) + ^m^ + (P) + {U{r)) + {Hbf) ■ (27) 

The static potential U{r) takes the place of L, the universal quark binding 
function, in Eq. 1 of Ref ^ . 



4 Mass Splitting Formulae 



The static potential in which the quarks interact determines the form of the 
NRQM wavefunction. The strong Coulombic term gives the largest energy con- 
tribution of terms within the Breit-Fermi interaction, and therefore would also be 
expected to substantially alter the wavefunction in perturbation theory. There- 
fore, we include the strong Coulombic term in the static potential: 



4 (T 

V{t) ^ U{r) - 

6 r 

Then the mass formula Eq. ^ becomes, using Eq. ^ 



(2^ 



M2 + {p) + Jm^ + + {V{r)) + aQqB 



+ aQq 



-a, 



2mM 2 



1 1 



■ (29) 



Now at last we are in a position to write explicit formulae for the mass 
splittings of interest. Denoting the mass of a meson of spin S and valence 
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quarks Q,q as M^{Qq), we define: 





= M\Qu) 


-M\Qd) 




= M\Qu) 


-M\Qd) 


A* 


= M\Qu) 


- M\Qu) 


A* 


= M\Qd) 


- M\Qd), 



(30) 



where u and d, the up and down constituent quarks, are nearly degenerate in 

<^ 1. Therefore, 



mass: defining Am = m„ — and m = 



Am. 



we have 

the differences in Eq. |30| are expanded in Taylor series in ^ about m. It is also 
convenient to define 



A = 



1 



fi = m(3 
Do.. ^ P 



1 + m/M' 
usual reduced mass, 



/i das 
as 9/i 



X = A,B,C,D,{V{r)) 



(31) 



Then the expressions for mass splitting are 



A 



0,1 



2m^ + Da 



+ 



Am ^ Am 
2m m 



— agAm 
3 



[Dc-C + CD^J 



IT 



2m^ 
+aQ 



, ^ m m 



(Dd + BB.,) - 2 1 + 25 



m 
M 



D 



B + 
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+0 



2mM 

Am-- 



:C 



m 



+ o \ a 



TT ^m m \ ^ 



Am 



m 



(32) 
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Note that no derivatives appear in the uem terms because we take both a em 
and (but not a^) as expansion parameters. Furthermore, the running of 
Q;s(/i) is exphcitly taken into account. 

For vector-pseudoscalar sphttings, we have 



A 



57r 



Qq 



3mM 

+0 



^ . 4 Am 
-as - aQq ) D ± 3"^-^ 



4 
3 

Am" ^' 



{Dd + DDo,^ - D) 



m 



+ o\a 



Am 



m 



^ , with ± for g = u{d). 



(33) 



Let us remind ourselves of the physical significance of the terms in the previous 
two equations. Terms containing A signify kinetic energy contributions, includ- 
ing intrinsic quark masses. The potential term is identified, of course, by V^; B, 
C, and D denote static Coulomb, Darwin, and hyperfine terms, respectively. 



5 Quantum-mechanical Theorems 

In order to apply the foregoing results, we will need to evaluate the expectation 
values A,B,C,D, and {V{r)) for our potential V{r). There are two quantum- 
mechanical theorems which make the evaluation of these expectation values and 



their mass derivatives simpler The first is the 



Theorem 1 (Feynman-Hellmann Theorem) For normalized eigenstates of 
a Hamiltonian depending on a parameter X, 

^^/^V (34) 

dx \ dx / ^ ^ 

In the particular case that X = fi, 

,35) 

Ojji jJ, I 

The other result may be less familiar. For reasons that will become clear, 
let us call it the 

Theorem 2 (Generalized Virial Theorem) Consider hound eigenstates 
ue{r) in a spherically symmetric potential V(r) such that 

lim r'^Vlr) = 0. 
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Then, writing the Schrddinger equation as 



in + 

and defining at by 
then 

i) ae is a nonzero constant; 
a) for q > -2i, 



E - V{r) - 



Ui{r) = 0, 



ut{r) _ 



lim 



2/i 



(2^ + l)XX-2/ = -~^lr'^-^\2q{E-V{r))-r^] 



dV\ 



dr I 



2£(£+l)_-5(5_2) 



(36) 



Clearly this theorem will prove most useful for potentials easily expressed 
as a sum of terms which are powers in r. But in fact there are some interesting 
general results included. For example, the q — £ — Q case generates the well- 
known result for s-waves, 

11 I dV \ 



i*(o)r 

whereas the q = 1 case produces 



27r?i2 \dr r 



dV 



(37) 



(38) 



the quantum-mechanical virial theorem. 

Using partial integration, the Schrodinger equation, and the GVT, it is 
possible to show the following {h — 1): 

A = 2pi{E-{V{r))), 



'duii/f 
dr 



C 



D 



dr 



4ti 



r I \ r 



1 /^^\ f ^' 



/i IdV 
27r \ , 



(39) 
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In addition, we must also uncover what we can about the /i-dependence of 
expectation values. For a general potential this is actually an unsolved problem. 
However, unless the potential has very special /i-dependence, it can be shown 
that only in the case V{r) = Vqt^ is it possible to scale away the dimensionful 
parameters Vq and yU in the Schrodinger equation. In that case, the //-dependence 
will be entirely contained in the scaling factors, and computing Dx will be trivial. 
Unfortunately, in the potential we consider in the next section, we will see that 
this is not the case, and we must resort to subterfuge to obtain the required 
information. 



6 Example: V{r) = ^ — - 



The potential V{r) = ^ — ^, where k = is interesting because it phe- 
nomenologically includes quark confinement via the linear term. This potential 
was considered in greatest detail by Eichten et al. [0 to describe the mass 
splitting structure of the charmonium system (and was later applied to bot- 
tomonium). The Schrodinger equation was solved numerically; currently, no 
analytic solution is known. However, it is possible to extract a great deal of 
information from their tabulated results, as we shall see below. 

This is possible because of the GVT. If we rescale the Schrodinger equation 
with the linear-plus-Coulomb potential to 

w 



^(^ + 1) A , \ , , 



0, 



where 



'2/i 



P 



1/3 



(40) 



then the GVT gives 



A = K{2iiaf'^ 
we{p) = Ui{r) 




1/6 



(41) 



(g = 0) al So^i 



2/i 



1 + A 



3 (p) - 2C 




(42) 
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Also, defining 



we find 



. dp 



{v') = -{p)+C + X 



.Pi 



j dp, 

iii + 1) 



(43) 



(44) 



It is a happy accident of this potential that all of the quantities in the 
expectation values we need, for any i, may be expressed in terms of the three 
quantities (, (^), and {v'^), which are exactly those values tabulated for the 

/2u\ 1/3 

Is-state, as functions of A, in Eichten et al. Table I. Defining a = {-^\ 
taking £ = (as per our mesonic model), we find 

A 



and 



2/2 

a [v 



B 
C 

D 



a 
2A 

Air 



3 v^)-C 



2BC + a(^-3 + \(^y^ 



X 



+ 1 



(45) 



So now we can compute all of the necessary expectation values numerically. The 
superficial singularity in B{\ = 0) is false; -B(O) is computed by extrapolation 
of the computed values of B for nonzero A and is found to be finite. 

The mass derivatives must be handled in a different fashion. We begin by 
defining 



n - 9C 



/i- 



and 



djji 

~ ^ /i da^ 

Us op 



d 



p- 



dp 



Prom the Feynman-Hellmann theorem (Eq. ^5]) we may show 



C 

3 



(46) 



(47) 
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As mentioned in the previous section, scaling of the Schrodinger equation can 
be accomphshed for //-independent potentials that are monomials. In the case 
A = (a purely linear potential), the scaling would be perfect, and (, (^) , 
and (f^) would be //-independent. In the A 7^ case, the derivatives must be 
found numerically. Again, fortunately, we have a table of numerical values of 
the desired expectation values, as a function of A(/i). We fit the expectation 
values Y (= (t-/ , {v'^)) to the functional form 



r(A) = Fo + 



Then, using Eq. ^ we find 

Dy 

Finally, define 



riY {Y - Fo) 



D 



X 



u--— for X 

on 



A,B,C,D, 



so that 

Then we find 
2 



13 D 



X 



(48) 



(49) 



(50) 



(51) 



Da 



D 



B 



-A + a'D^, 
2A 2 



(A^O), 



D 



c 



D 



D 



\c + 2(7^ I (-C + ^c) ^ + C^B + <y 



- + D^\l-)+D, 





(52) 



In the exceptional case of Db-, we simply note that, for A = 0, we have perfect 
scaling of the wave equation, and we can quickly show that Db ^ ^ = ^ ^. 
This provides us with everything we need to produce numerical results. 

Before leaving the topic, let us mention that many complications of /i- 
derivatives of expectation values vanish if the potential itself has the appropriate 
/{-dependence, for then scaling of the wave equation is possible. For example, 
one can scale the Schrodinger equation for the potential 



V{r) 



cjjL r , 

r 



(53) 
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where c is a pure number. 



7 Numerical Results 

The method of obtaining results from the theory requires us to choose several 
numerical inputs, most of which are believed known to within a few percent. 
Let us choose the following inputs to the model: 

m = 340 MeV, M, = 540 MeV, 
Mc = 1850 MeV, Mb = 5200 MeV, a = 1.95 GeV"^ (54) 

The light quark constituent mass is arrived at by assuming that nucleons consist 
of quarks with negligible anomalous magnetic moments, which can be added 
nonrelativistically to provide the full nucleonic magnetic moment. Likewise, 
the strange quark mass issues from the same considerations applied to strange 
baryons . The c and b quark masses are simply found by dividing the threshold 
energy value for charm and bottom mesons by two (however, smaller masses have 
been predicted using semileptonic decay results in addition to meson masses 
||13||). The confinement constant is inferred from charmonium levels |ll2[ . 

One important input not yet mentioned is Am, the up-down quark mass 
difference. Traditionally, this assumes a value of ~ —3 to —8 MeV, in order to 
account for the electromagnetic mass splittings of the lighter hadrons. In this 
model, with the inputs listed in Eq. ^ we find that the experimental splittings 
for the D- and i?-mesons (both vector and pseudoscalar) can be satisfied within 
one standard deviation of experimental error for values of Am in the narrow 
range of —4.05 to —4.10 MeV. In contrast, it is found that for no choice of Am 
can one simultaneously fit D- and ii'-meson data simultaneously, as was done 
in the earlier models. 

Before exhibiting the quantitative results, let us describe the method by 
which they are obtained. Once particular inputs for the above variables are 
chosen, one can compute the various mass splittings for the values of A oc that 
occur in Table I of Ref. []T2| , and in-between values may be interpolated. We then 
fit vector-pseudoscalar splittings (computed via Eq. to the corresponding 
experimental data (since these numbers have the smallest relative errors of the 
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splittings we consider) and thus obtain a value of a^. For the three systems K, 
D, and B, we use the three values of as to estimate graphically (and admittedly 
rather crudely) its mass derivative. Applying the values of the strong coupling 
constant and its derivative to the splittings in Eq. we generate all of the 
other values. If the resultant numbers do not fall within the experimental error 
bars for such splittings, we vary the input parameters (most importantly, Am) 
until a simultaneous fit is achieved. 

Table I displays the various contributions to mass splittings derived in this 
fashion for B- and Z^-mesons. Although the kinetic term (which includes the 
explicit difference Am) and the static Coulomb term are unsurprisingly large, a 
significant contribution to the mass splitting arises in the strong hyperfine term. 
That strong contributions to the so-called electromagnetic mass splittings could 
be important was observed by Chan 0, and was exploited in the subsequent 
literature. It is exactly this term which is most significant in driving the B 
splittings toward zero. Note also the decrease in the derived value of as as the 
reduced mass of the system increases when we move from the D system to the B 
system, consistent with asymptotic freedom in QCD. It was this running which 
motivated the inclusion of mass derivatives of the strong coupling constant in 
this model. If they are not included, one actually obtains a value of Am > 0, in 
contrast with all estimates from both nonrelativistic and chiral models. 

The net result is that one can satisfactorily fit the data for the D and B 
systems simultaneously in the most natural nonrelativistic model with a physi- 
cally reasonable potential. The comparison of the results of this calculation for 
Am = —4.10 MeV to experimental data is presented in Table II. 

However, the table also exhibits very bad agreement for the K system 
(despite the fact that the fit to vector-pseudoscalar splittings yields the value 
as = 0.424, which runs in the correct direction). One may view this as a failure of 
the nonrelativistic assumptions of the model in a variety of ways. Most obvious 
are the ansatz Eq. which is certainly not an airtight assumption in even the 
best of circumstances, and the crudeness of the estimate of Other possible 
problems include the assumption that the quarks occur only in a relative i = 
state (relevant for ii'*-mesons) , and the assumption that the strong effects are 
dominated by a confining potential and one-gluon exchange, since at the lower 
energies associated with the K system, o(a^) terms and more complex models of 
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confinement may be required. The failure of tliese assumptions can drastically 
alter the strong hyperfine interaction, which determines the size of ag-, and hence 
the other mass splittings. 

Some may find the small size of somewhat puzzling. This is primarily the 
result of the confining term of the model potential: it causes the wavefunction 
to be large at the origin, and thus a small is required to give the same 
experimentally measured vector-pseudoscalar splitting (see Eq. |33| ). Such small 
values for the strong coupling constant might lead to excessively small values of 
^QCD and large values for mesonic decay constants /g^. Indeed, given the naive 
expressions for these quantities: 

^ (33-2.iTlog(«l)' 

and, assuming the relative momenta of the quarks is small, 

1 2 

let us consider, for example, the D system. Then = 0.363 and /i = 287 MeV, 
and with three flavors of quark, we calculate Aqcd = 42 MeV and fs = 342 
MeV. However, one may state the following objections: first, Aqcd is computed 
from the full theory of QCD, but the nonrelativistic potential approach includes 
the confinement in an ad hoc fashion, by including a confinement constant a, 
which is independent of a^. Furthermore, choosing Aqcd as the renormalization 
point forces an artificial singularity at /i = Aqcd] the problem is that little is 
known about the low-energy behavior of strong interactions. At low energies the 
computation and interpretation of Aqcd requires a more careful consideration 
of confinement. With respect to the decay constant, the assumption that the 
quarks are relatively at rest leads to the evaluation of the wavefunction at zero 
separation. Inclusion of nonzero relative momentum will presumably result in 
the necessity of considering separations of up to a Compton wavelength r ^ j^, 
for which the wavefunction is smaller in the Is-state. Thus decay constants may 
be smaller than computed in the naive model. 

There is one further qualitative success of this model, a partial explanation 
of the experimental facts that D* - D, = 141.5 ±1.9 MeV Ri D* - D, and 
B* — Bg = 47.0 ± 2.6 MeV B* — B, namely, the approximate independence of 
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vector-pseudoscalar splitting on the light quark mass. In our model, the leading 
term of the splitting is, using Eqs. |33| and H3, 
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9Ma2 



A 



+ 1 



(57) 



Inasmuch as /3, A {^^y, and as are slowly varying in the light quark mass m, the 
full expression reflects this insensitivity, in accord with experiment. In fact, we 
may fit the experimental values above to obtain more running values of a^: 



141.5 MeV for a, = 0.351, 



A^, = 47.0 MeV for a, = 0.295, 



(5^ 



and again these decrease as the mass scale increases. Note, however, one kink 
in this interpretation: the heavy-strange mesons all have larger reduced masses 
than their unfiavored counterparts, yet the corresponding values of are nearly 
the same. 



8 Conclusions 

In this paper, we have seen how mass contributions to a bound system of par- 
ticles are derived from an interaction Hamiltonian in field theory, and how this 
calculation is then reduced to a problem in nonrelativistic quantum mechanics. 
For the system of a quark and antiquark bound in a meson, the exchange of one 
mediating vector boson reduces to the Breit-Fermi interaction in the nonrela- 
tivistic limit. It is also important to consider contributions to the total energy 
from the kinetic energy and the long-range potential of the system; in fact, the 
higher-order momentum expectation values can be so large that it is necessary 
to impose an ansatz (Eq. |26|) in order to estimate their combined effect. Future 
work may suggest better estimates. 

It is found in the case of a linear-plus- Coulomb potential that the largest 
contributions to electromagnetic mass splittings originate in the kinetic energy, 
static Coulomb, and strong hyperfine terms. However, it is likely that simi- 
lar results hold for other ansatze and potentials. As in other models, vector- 
pseudoscalar mass differences are determined by strong hyperfine terms. 
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With typical values for quark masses, the confinement constant, and the 
up-down quark mass difference, we can obtain agreement for the mass split- 
tings of the D- and S-mesons. The failure of the model for K mass splittings 
is attributed to the collapse of the nonrelativistic assumptions in that case. 
The model also qualitatively explains the similarity of heavy-strange to heavy- 
unfiavored vector-pseudoscalar splittings, although additional work is needed to 
explain why these numbers are nearly equal, despite the expected inequality of 
as at the two different energy scales. 

Another interesting problem is the running of ctg itself at low energies. As 
mentioned in the results section, this running cannot be neglected if we are to 
obtain sensible results, and yet our approximation of this running is based on 
crude assumptions. The size of ag also enters into another possible development, 
namely, whether terms of o(q;^) are important, particularly for the K system. 
More reliable estimates are required. 

In addition to the explicit formulae derived in this paper, the techniques 
employed here may be applied to later efforts: in particular, the explicit consid- 
eration of the mass-dependence of expectation values and the use of quantum- 
mechanical theorems to relate various expectation values for certain potentials. 
The methods and formulae in this work may prove to be a starting point for 
subsequent research. 
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Table I: Contributions to mass splittings of heavy m esons 





D mesons 


B mesons 




0.363 


0.312 


Source 


MeV 


MeV 


Tt;nt;rii n r»j^ ir^ 






Kinetic energy 


-4.109 


-3.523 


Potential enersrv 


1.057 


-1.645 


Strong Darwin 


-0.834 


-0.635 


EM Darwin 


-0.769 


0.147 


Stritip dmilomli 


-2.442 


1.252 


Q 






Strong hyperfine 


2.148 


4.075 


EM hyperfine 


0.424 


-0.561 


Total 


-4.525 


-0.889 








Strong hyperfine 


3.683 


5.244 


EM hyperfine 


1.817 


-0.825 


Total 


-1.596 


0.017 


1- -0- 






Strong hyperfine 






(leading) 


141.30 


46.04 


(subleading) 


± 0.77 


± 0.58 


A* 






EM hyperfine 


0.93 


-0.18 


Total 


143.00 


46.45 


A* 






EM hyperfine 


-0.46 


0.09 


Total 


140.07 


45.54 
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Table II: Meson mass splittings compared to experiment 



Mass splitting 


Notation 


Predicted (MeV) 


Expt. (MeV) 




A° 


-0.98 


-4.024 ± 0.032 






-0.15 


-4.51 ± 0.37 " 


K*+ - K+ 


A* 

su 


398.6 


397.94 ± 0.24 " 


K*^ - 




397.8 


398.43 ± 0.28 " 




A° 


-4.53 


-4.77 ± 0.27 




A^ 


-1.60 


-2.9 ± 1.3 


D*^ - 


A* 

cu 


143.0 


142.5 ± 1.3 


D*+ - D+ 




140.1 


140.6 ± 1.9 " 




K 


-0.89 


-0.1 ± 0.8 


^*+ — 




0.02 


NA 


B*+ - B+ 


AL 


46.5 


46.0 ± 0.6 ^ 


B*^ - B^ 




45.5 


46.0 ± 0.6 ^ 



" obtained as a difference of world averages 
average of charged and neutral states 
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Figure Captions 

FIG. 1. Diagrammatical representation of Al/j. 

FIG. 2. Diagram for M.fi in the mesonic system. 

FIG. 3. Free quark Feynman amplitude Ai. 

FIG. 4. Notation and conventions for the mesonic system. 
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Pf 



Q, helicity H 



q, helicity h 



Pi 



Pi^^-P 



2 
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